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SUMMARY 


The flow field of three-dimensional incompressible wall jets 
prototypic of thrust augmenting ejectors with large cross flow is 
solved using a very efficient centered-Euler scheme in an orthogonal 
curvilinear coordinate system. The computational model treats initial 
conditions with arbitrary velocity profiles at the jet exit. An 
averaging approach is employed for the first few marching steps to 
overcome spurious numerical oscillations associated with arbitrary 
initial profiles. Laminar as well as turbulent wall jets are simulated. 
Turbulence is introduced using a two layer mixing length model appropriate 
to curved three-dimensional wall jets. Typical results quantifying 
jet spreading, jet growth, nominal separation and jet shrink effects 
due to cross flow are presented. 


*This work was sponsored by the Naval Air Development Center under 
Contract N62269-77-C-0412 . The monitor for this effort was 
Dr. Kenneth Green. 



INTRODUCTION 


Modern naval aircraft can reduce strike force vulnerability by 
the attainment of vertical lift-off capability. To achieve accelerations 
associated with typical payloads, a high augmentation ratio 4) is 
required. Various propulsive lift concepts have been advanced toward 
obtaining this goal. In the XFV-12A, an ejector system composed of 
a centerbody and two Coanda wall jet-s is currently under development. 

A central feature of the flow fields produced by this device is three 
dimensionality. This has been particularly evident in subscale flow 
visualization on the Coanda surfaces. It is believed that these flow 
processes may be important toward 4i maximization. One way of understanding 
this relationship is through theoretical modeling which can provide a 
means of reducing the high cost of powered lift testing. Unfortunately, 
existing methodology has been limited in the past to two-dimensional 
flows for the analysis of wall jets and complete ejector systems. 

Analytical methods and computational algorithms are therefore necessary 
to compute three-dimensional flows typical of reality. 

To shed light on typical flow patterns encountered, due to the 
effect of taper and sweep on augmenter wings as well as upper-surface- 
blown configurations, a study, "Three-Dimensional Flow of a Wall Jet," 
was initiated by the Naval Air Development Center to investigate wall 
jet flows which exemplify typical features of more complex propulsive 
lift applications. The purpose of this study has been to apply modern 
computational methods to the treatment of wall jet flows with three 
dimensionality . 

The formulation employs boundary-layer equations in an orthogonal 
curvilinear coordinate system. It can be shown from a systematic 
order of magnitude analysis that the boundary-layer equations also 
apply for wall jets, providing the distance from the jet exit is 
sufficiently large to establish complete mixing, the Jet height is small 
compared to a characteristic radius of curvature, and the Reynolds 
number based on the exit height is large. A transformation is Incorporated 
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to stretch the coordinate normal to the flow. At streamwise planes, 
the resulting nonlinear partial differential equations are treated as 
ordinary differential equations. These are solved using a very efficient 
two-point boundary value finite-difference method devised by Keller 
and Cebecl^”^ known as **box method.” The turbulence is introduced using 
a two-layer mixing length model appropriate to three-dimensional wall 
jets. 


Equations in Curvilinear System 


The governing equations for three-dimensional incompressible flows 
over a wall jet in a curvilinear orthogonal coordinate system shown in 
Figure 1 are given by the following equations: 


Continuity 


0/3x)(h^u) 4 - 0/3z)(h^w) + 0/3y)(h^h^v) = 0 


( 1 ) 


x-Moraentum 
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z- Momentum 


u 9w 9w ^ ^w 


uwK + u^K 
2 1 


+ A 

ph^ 9z 9y 




(3) 


Here h^ and h^ are metric coefficients and are functions of x and z, 
and the parameters and K 2 are known as the geodesic curvatures of 
the curves z = const- and x =■ const., respectively. 
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Figure 1.- Physical system and flow schematic. 
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The boundary conditions for Eqs. (1) through (3) for zero mass 
transfer are 


y = 0 


u,w,v = 0 


y oo U u^(x,z) W ^ w^(x,z) (4) 

As indicated earlier, the previous equations are transformed by 
defining 

x*x z = z r)» (u^/vs^)V 


and introducing a two-component vector potential given by 


h^u 


3y 


h^w 


3y 


^iV 



( 6 ) 


In addition, the dimensionless variables f and g related to and (}) are 
defined by 



(7a) 


't> = 


(u vs 
e 




e e 


(7b) 


Here s^, which denotes the arc length along the x coordinate, is 
defined by 


s 


1 



( 8 ) 
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The parameters and in Eq. (2) are defined by 


1 

If = _ ^ L 

’ dz 


and 


^2 = 




dh 

^ 

dx 


With the concept of eddy viscosity and with the previous transformed 
variables, it can be shown that the system of Eqs. (1) through (4) can 
be written as 


x~Momentum Equation 

(bf")’ + Pjff" + ?^[l - (f’)2] + p^[l _ f.g.] 
+ Pgf’g + P3[l - (g’)2] 


ft ML 

- f" + p / 

g’ -- 

f.t i&V 

9x 

9x 7 \ 

* 3z 

9z / 


z-Momentum Equation 
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n = o f*g"f'“g'=o 
n = f’ = g' = 1 


(9) 


( 10 ) 

(lla) 

(llb) 
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Here the primes denote differentiation with respect to 


and 


b = l+e"*" *■ e/v f = g' = w/w^ 

The coefficients P, to P, . are functions of u , w , h, , h, , K , and 
110 

and are given by the following formulas : 


Pj = (Mfl)/2 - SjK^ ?2 = M P 3 = R 
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(12a) 


(12b) 


In order to solve Eqs. (9) through (11), initial conditions are 
required at a starting plane. In the case of the boundary-layer problem. 



the initial conditions at x = 0 and z = 0 planes are obtained by solving 
the limiting form of Eqs. (9) and (10). For a wall jet, initial velocity 
profiles are prescribed at some downstream x plane and along the z = 0 
plane, attachment line equations are solved. The attachment line equations 
are obtained by differentiating the z-momentum equation with respect to 
z and setting 


3p ^ 3u ^ 3v ^ 

3z 3z 3z 3r^ 


0 


The resulting attachment line equations valid at the z = 0 plane are 

(bf’)’ + p ff" + p ri - (f’)M + p,gf" = xp, f 

1 2 3 ^ 1 0 I 

(bg")' + p^fg" + p^(i - f'g') + pji - (g')2] 



... ' 

3x 

■ ^ 3l_ 


(13) 


+ P^gg" - xPj, 


f - l£l _ li 

9x ® 3x 


(14) 


Here, g' is defined as w /w 

z e 


Eddy-Viscosity Model 


Eqs. (2) and (3) contain Re3molds shear stress terms -u'v’ and -v’w*. 
In order to satisfy the closure assumptions for these shear stress terms, 
we use the eddy viscosity concept and define 


1 - 


1 + <2 y “y 


1 - 


C W < T 
2 1 _ 

1 + K y w 
1 ^ y 


(15) 
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The second term inside the bracket in Eq. (15) is due to curvature where 
K, and denote the radius of curvature of z = const, and x = const, lines. 

i 2 

The quantity e is assumed to be same in both the x and z directions and 
is represented by a two-layer model. Referring to Fig. 2, the structure 
of these layers is as follows: 


First Layer 

e = (0.435 y)^ ^ u^ + 0 < y ^ y* 

Second Layer 

e = (0.125 y ^u^ + w^ y > y* 

l y y 

where at y =» y^. 




j 


I - ^u^ + w^ I 

® ^ ■ ■ ■ ^ 0.01 

v/u' + w' 
e e 


(16) 


and y* is obtained by imposing continuity in e at y = y*. This yields 
y* = y^ . and appearing in Eq. (15) can be assigned values 
between one and three. 


Finite Difference Equations 

First, reduce the system (9)-(10) to the first order system 


f = u 

(17) 

u' = V 

(18) 

g' = w 

(19) 

w' = t 

(20) 
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Figure 2.- Two-layer eddy viscosity model. 
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(bv)' + Pjfv + Pj(l - u^) + P,[l - u»l 


+ P,vg + P,[l - w'l - xP^, 


f 9u 3f / 9u ^ 

^ ^ 9l^ ^ ^ 


( 21 ) 
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14 3 6 
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Then, using the box method, we have 
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(24) 
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We use the notation 



and 


P 


n-h 

i-h 


V , 

J 


u . 
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Equation (21) becomes, 


with the box centered at 




) 
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(27) 


Equation (22) and the attachment line equations (13)~ (14) are discretized 
similarly. Details of the procedure are given in Ref. 4. 


152 



The solution procedures Involves the following steps: 

(1) Solve the attachment line equations (13)-(14) with boundary 
conditions (11) at x = and z = 0 assuming initial conditions 
on X = Xq . 

(2) March in the z-direction along the plane x = x^ and solve 
equations (17)-(22) with boundary conditions (11) for the 
unknowns (f,u,v,g,w,t). 

(3) Repeat steps (1) and (2) for the next x-plane, x = x^ , and so on. 

The most efficient way to solve the finite difference equations is 
to use a pseudo-Newton's relaxation scheme. These equations may be written 
as a system of nonlinear algebraic equations by writing 


$(u) =* 0 


where 


u = (f 


n,i n,i n,i n,i n,i 


w. 




1-0 


Then, the relaxed Newton's method becomes 


3u 


5u 


(v-1) 



(28a) 


u 


(V) 


(v-1) 


+ (jj6u 


(v-1) 


(28b) 


for V = 1,2,... 
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The method is said to have converged when 


<e 

^ oo 


(a prescribed error tolerance) 


We call Eq. (28) a pseudo-Newton's method because we linearize the 
b terms in equations (21) and (28) by evaluating them at v— 2 before 
computing the Jacobian matrix, 3$/3u. Consequently, this algorithm will 
not be quite quadratically convergent. We, therefore, employ relaxation 
^ 1) to accelerate it. Remarkably, underrelaxation (u) < 1) works 
very well, while overrelaxation (oo > 1) diverges. Values of w of 0.5, 

0.6, 0.7, 0.9, and 0.9 all give good results with oo = 0.7 slightly the 
overall best for some of our computational experiments. 

An important feature of Keller's box method is that the Jacobian 
matrix can be put into block tridiagonal form and very efficient elimination 
schemes can be employed for solving equation (28a) . 


Minor Difficulties with the Numerical Algorithm 


When starting at x =• x^ 0 with supplied velocity profiles, 

unnatural oscillations developed in the solution. This difficulty was 

eliminated completely by employing the following "trick." The first 

10 mesh points in the x-direction were set at k = 10~‘*. For the first 

n 

five planes in the x-direction and all points in the z-direction in these 
planes, an average value was used for past points, i.e.. 




0.5 
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Beginning with the sixth x-plane, the averaging was eliminated (the standard 
algorithm was employed) . At the eleventh x-plane a geometric mesh- 
stretching algorithm of the following form was used: 

k - 1.2k , , n = 11,12,13,... 

n n-1 

No such stretching has been employed in the z-direction, but in the future 
it may also be required for rapidly changing profiles. It should be 
noted that our averaging algorithm was required in both the x and z 
directions to remove all oscillations. 

A mesh refinement algorithm is used which adds or deletes points 
depending on the relative local variation in the truncation error of the 
difference equations. Roughly 80 grid points in the f|— direction and 
11 grid points in the z-direction are employed. 

Results 


Computations were performed on the Berkeley CDC 7600 machine. A 

typical calculation required about 6 minutes of CPU time. Fig. 3 indicates 

the external and Initial velocity distributions which have been used as a 

basis for our calculations. The parameter 0 was introduced as shown to 

vary the initial cross flow while keeping the total velocity constant 

as a rough simulation of a fixed supply of engine mass flow. The velocity 

profile was selected to have a characteristic fully developed character 

associated with turbulent wall j^t flows. Future aspects of this effort 

will consider the "eating up" of the potential core which is assumed to 

occur upstream of the Initial station of this analysis. The parameters 

C and C were chosen to provide slope and value continuity of the profile 
2 . 3 

at V = V For v ^ V the profile has a half Gaussian character 

^ ^ max max 

associated with a free jet. For y < y the profile has a boundary- 

IIL3X 

layer character. In the examples, the u and w initial profiles were 
assumed to be identical. Moreover, the 6 distribution was selected to be 
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EXTERNAL VELOCITY 
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Figure 3.- Initial and boundary conditions. 



qualitatively similar to that observed by rake surveys on the XFV-12A. 

The zero cross flow case was achieved by setting C to 10 

Figure 4 demonstrates decay of the peak velocity with the standardized 
distributions of Fig. 3, with and without cross flow. It is evident from 
the figure that cross flow has a dramatic effect on enhancing the decay 
of the maximum velocity. In the calculations, the exponent n in the 
external velocities is assumed as Hy roughly in accord with a value 
obtained from a two-dimensional line sink simulating inflow originally 
proposed by G.I. Taylor.^ In a more realistic model, these external 
velocity distributions should be corrected for three-dimensionality and 
elliptic interaction with the wall jet. A calculation of this type 
would be a more accurate representation than the present approach of 
planform and surface curvature effects. In this connection, we recognize 
that that means of simulating taper, sweepback, and spanwise pressure 
gradients in the present analysis is solely through cross flow adjustment. 

The three-dimensional inviscid potential (p can be characterized by 
a surface sink distribution of the form (see Fig. 5) 


1 /7* q(g>C)dgdC (29) 

-^-4 y (x-C)^ + + (z-O 


where S the area of integration refers to the total jet area on and 
off the wing. The quantity o is the sink strength obtained by matching 
with an "outer limit" of the second order solution for the velocity 
normal to the body appearing in the viscous inner wall jet solution. 

The quantity o for two-dimensional boundary layers is analogous to the 
streamwise gradient of the displacement thickness 6'(x). To Include 
]_iftlng surface effects, a surface doublet or vortex distribution should 
be added to (29). The local vortex strength can also be determined by 
matching. 

The inflow velocity related to the sink intensity O in (29) is in 
turn a function of the entrainment. This quantity is also significant 
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Figure 5.- Tapered thrust augmented wing (TAW). 
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from the standpoint of the tradeoff between skin friction, BLC, and rapid 
acceleration of the secondary in compact three-dimensional thrust augmenting 
ejectors such as those employed on the XFV-12A and upper surface blowing. 

In Fig. 6, the comparison between cross flow and the absence of it 
gives the indicated entrainment variations with streamwise distance. In 
spite of the appreciable Increase in decay of the maximum value of u shown 
in Fig. 4, and resultant shear stress in Fig. 7, only a slight difference 
in entrainment quantity and rate is shown in Fig, 6. The difference in 
maximum velocities which are similar for w, the spanwlse component, are 
presumably related to the enhanced dissipation associated with cross flow 
and that implied by the eddy viscosity model. The lack of a corresponding 
decrease in entrainment rate may be due to nonlinear compensating effects 
built into the turbulence model and cannot be readily explained on an 
intuitive basis at this time. In this connection, other calculations 
will be performed for which the streamwise component of the initial 
velocity is held fixed rather than its overall magnitude on introduction 
of cross flow. It should be noted that the expression for entrainment Q 
given in Fig. 6 assumes that w at the tip z = ^tip ” this is not 

the case, an additional term must be added to this relation. 

Associated with the previous results. Fig. 8 shows the effect of 
cross flow on jet spreading rate related to previously, only 

small differences are indicated for the cases considered. In Fig. 9, 
however, an important upstream movement of the separation line is indicated 
with the introduction in cross flow. This result is significant with 
respect to penalties associated with taper and sweep in three-dimensional 
ejector diffusors. 

In Fig. 10, another important consequence of cross flow is examined 
in connection with the surface streamline pattern. In the figure, two 
cases are compared involving differing amounts of cross flow. Significant 
enhancement in downstream streamtube contraction is obvious with increase 
in cross flow. This contraction could presumably lead to end wall 
separation of the type observed on the XFV-12A. 
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Figure 10.- Cross flow effect on jet ^shrink" and ’’end-wall pull away 




Conclusions 


A class of cases were investigated roughly possessing Initial flow 
angularity and adverse pressure gradients prototypic of those on the 
XFV-12A. Results obtained from the computational model indicate that if 
the initial total velocity is kept fixed then the introduction of the cross 
flow enhances the decay rate of the peak of the streamwise velocity 
component. In addition, the entrainment quantity and its rate decrease 
with increased cross flow. The implication of this phenomenon with 
respect to taper effect on boundary layer control (BLC) of the XFV-12A 
Coandas is not as significant as a "jet shrink" which has also been 
indicated in our approximate three-dimensional model. This contraction 
has been postulated as a mechanism promoting end-wall separation. To 
our knowledge, our model is the first to quantify such trends. Finally, 
the effect on the prescribed external adverse pressure gradient in the 
presence and absence of cross flow has also been examined. From the 
limited results, the spanwise separation line moves progressively further 
upstream with increasing cross flow. 
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